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Abstract. A particular continuous-time multitype branching process is consid- 
ered, it is the continuous-time embedding of a discrete-time process which is very 
popular in theoretical computer science: the m-ary search tree (m is an integer). 
There is a well-known phase transition: when m < 26, the asymptotic behavior of 
the process is Gaussian, but for m > 27 it is no more Gaussian and a limit W of 
a complex- valued martingale arises. Thanks to the branching property it appears 
as a solution of a smoothing equation of the type Z — e~ XT (Z^ + ■ ■ • + Z^), 
where A e C, the are independent copies of Z and T is a R+- valued random 
variable, independent of the Z^ k > . This distributional equation is extensively stud- 
ied by various approaches. The existence and unicity of solution of the equation 
are proved by contraction methods. The fact that the distribution of W is abso- 
lutely continuous and that its support is the whole complex plane is shown via 
Fourier analysis. Finally, the existence of exponential moments of W is obtained 
by considering W as the limit of a complex Mandelbrot cascade. 
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1 Introduction 

Consider a continuous-time multitype branching process (X(t),t > 0). Types are 
seen as colors of particles and there are m — 1 colors, where m > 2 is an integer. 
The reproduction of the process is given by a particular matrix R (written in 
( 12. II) ). and any particle of colour j lives a random time of exponential distribution 
with parameter j. Such a classical process is considered for example in Athreya 
and Ney [1] or Janson [12] and it is precisely defined in Section 2. 

When it is stopped at the n-th jump time, this process is nothing but the 
composition vector process {X^ T , n > 0) say, of an m-ary search tree, which is 
an important algorithmic structure in computer science. A numerous literature is 
devoted to the asymptotic behavior of this composition vector. A famous phase 
transition appears. When m < 26, the random vector admits a central limit 
theorem with convergence in distribution to a Gaussian vector: see Mahmoud 
and Pittel [18], Lew and Mahmoud [13]. When m > 27, it has been proved that 

X° T = nv\ + ^(n X2 W DT v 2 ) + o(n a2 ) a.s., 

where A2 is a complex number having a real part in ]^, 1[ (it is an eigenvalue of 
the replacement matrix R), where i>i,i>2 are deterministic vectors and W DT is 



the limit of a complex-valued martingale. Heated conjectures about the random 
variable W DT remain open (see jl], [2D], Chern and Hwang [6], Mahmoud [17], 
Janson [T2]). 

This article is focused on the asymptotic behavior of the continuous-time pro- 
cess (X(t)). The links between continuous-time and discrete-time processes are 
detailed in Section 4. In particular, we give the so-called martingale connection 
that relates the almost sure limits of both processes. 

Inspired by the methods used for a two-color Polya urn in [5], we first prove 
in Section 3 that for m > 27, X(t) admits the following asymptotic expansion: 

X{t) = e'^i(l + o(l)) + $(e X2t Wv 2 )(l + o(l)) a.s. and in L p Vp > 1, 

where £ is a Gamma distributed random variable and W a C-valued one. 

We are interested in the limit random variables Wk, k = 1, . . . , m — 1, each 
corresponding to Xk(t) which denotes the process X{t) when it starts from one 
particle of color k. Using the branching property, a system of dislocation equa- 
tions is written for the random vectors Xk(t) in Section [57X1 A system of fixed 
point equations satisfied by the corresponding limit laws is then derived in Sec- 
tion 15.21 In particular, the complex-valued random variable W\ is a solution of 
the fixed point equation 

Z = e - A ^(z« + --- + Z (m) ), (1.1) 

where {Z^ : k > 1} are independent copies of Z, T — + ■ ■ ■ + T( m _i), 
{t(j) '■ j > 1} are random variables independent of each other and independent 
of {Z^}, each t^) has distribution £xp(j) (we denote by £xp(j) the exponential 
distribution of parameter j: it has density x H- je - ^ on ]0, oo[). 

Further properties of W\ are derived from a fine study of Eq. (11.11) . We first 
show in Theorems 16. 121 and 16. 141 that Eq. ( II. ip admits a unique square- integrable 
solution having a given mean. In particular, this implies that Eq. ( II. ip charac- 
terizes the distribution of W\. This result is proven by two contraction methods 
applied to the corresponding smoothing transformations. The first one deals with 
suitable spaces of probability measures where the classical Wasserstein metric is 
adapted to the complex field; it leads to Theorem 16.121 The second contrac- 
tion method, that gives a proof for Theorem I6.14| consists in working on Fourier 
transforms of solutions and provides a somehow simpler proof. Furthermore, 
this second method gives a result of existence and unicity for solutions of the 
convolution equation 

r+oo 

$(t) = / $ m (t - u )f T (u)du, t e C, 
Jo 
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in a convenient space of functions, where fa denotes the density of T (see Re- 
mark 16.161) . 

Once the characterization of W\ by Eq. ( II. ip is proven, it suffices to derive 
properties of solutions of this distributional equation. We show in this way the 
following results on the law of W\. 

Theorem 1.1 When m > 27, the complex-valued random variable W\ admits a 
density and its support is the whole complex plane. Its Fourier transform satisfies 

Ee i(tm = o(\t\- a ) 

when \t\ — > +oo, for some a > 1. 

This theorem is a direct consequence of Theorem 17.171 that provides such 
properties for solutions of fll.ip admitting a nonzero mean. Our proof consists in 
showing successively that the characteristic function of any solution has modulus 
equal to 1 only at the origin, that it tends to zero at infinity, and finally that it is 
of order 0(|t| _a ) as \t\ — > oo for some a > 1 so that it is square-integrable on C. 
In the approach we need to prove a non-lattice property of Eq. f ll.ip via Gelfand- 
Schneider theorem, using the algebraicity of A2 (see proof of Lemma l7.20p . 

Theorem 1.2 When m > 27, the random variable W\ admits exponential mo- 
ments in a neighbourhood of the origin of the complex plane. If L\(z) = Ke zWl 
denotes its Laplace series, then L\ is holomorphic near and, after a change of 
variable, the function z H- —-L\ (z~ X2 ) is a solution of the differential equation 

y(m-l) _ m 

Theorem 11.21 is immediately derived from Theorems 18.271 and 18.301 just as 
Theorem 11.11 was derived from Theorem 17.171 To prove Theorems 18.271 and 18.301 
we consider a solution of (11.11) as the limit of a complex Mandelbrot cascade. The 
results are a consequence of fine analytical properties of the Fourier transform of 
the limit variable. 

The paper is organized as follows. 

The continuous-time multitype branching process is defined in Section [2j Its 
relation with the m-ary search tree is detailed in Section [3], while Section [4] is 
devoted to the second order asymptotic expansion of (X(t)) t >o and to its connec- 
tion with the corresponding discrete process. In Section [5J we use the branching 
property of the process to show that the martingale limits of the continuous- 
time process are related by a system of equations in law so that the fixed point 
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equation (II .ip emerges. These first four sections constitute the first part of the 
paper. 

The second part of the paper consists in putting the focus on Eq. (11 .ip that 
turns out to characterize the distribution of W\ so that all results on solutions 
provide results on W\. In Section [6] we define the natural smoothing transform 
associated with Eq. (II .ip and we show that it defines a contraction in the space of 
square-integrable probability measures with given mean. Results on the support 
and on absolute continuity of solutions are obtained in Section [7J Finally, Sec- 
tion [8] is devoted to the exponential moments and the Laplace series of solutions. 

2 Definition of the branching process 

In this section we introduce the definition of the continuous time multitype 
branching process (X(t)), and present the spectral decomposition of its tran- 
sition matrix. 



2.1 Infinitesimal generator 

In the whole paper, the underlying vector space is M m_1 or sometimes C m ~ l . Let 
R be the following square matrix of order m — 1: 

/ - 1 1 

-1 1 



R 



■■. 1 

\ m — 1 / 



(2.1) 



and for k — 1, . . . , m — 1, let be the k-th row vector of R: when 1 < k < m — 2, 
the k-th coordinate of w k equals —1, the {k + l)-th equals 1 and all the others 
are 0; u> m -i has m as first coordinate, —1 as last one, and for all others. 

Let G be the operator defined on functions / from C" 1 " 1 to any real or complex 
vector space by the following formula: for any vector v in C m_1 , 



m—1 



G(f)(v) = klk(v)\f(v + w k ) - /(«)], (2.2) 



k=l 



where 4 are the coordinate forms: h(xi, . . . , x m _i) = 
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Definition 2.3 The right- continuous process X = (X(t),t > 0) is the only con- 
tinuous time Markov process with state space M m_1 having G as infinitesimal 
generator. 



Equivalently, X is a continuous time multitype branching process with m — 1 
types (or colors), having R as reproduction matrix. The fc-th coordinate of the 
vector X(t), namely lk(X(t)), is the number of particles of color k at time t. A 
particle of color k lives a random exponential time with parameter k; when it 
dies, it reproduces one particle of color k + 1 if k — 1, . . . , m — 2, and m particles 
of color 1 if k = m — 1. 

This branching continuous time process can be thought as the embedded 
process of a discrete Markov chain X DT = (X^ T ) n£ -M which is a Polya-type 
discrete Markov chain associated with the node process of an m-ary search tree, an 
important algorithmic structure in computer science. This connection is detailed 
in Section |3J 



2.2 Spectral decomposition 



Let Rq be the matrix of G's restriction to linear forms in the canonical basis 
{h)i<k<m-i- One immediately checks that 



/ 



Rg 



1 

-2 



\ 



2 

-3 



m — 2 
— (m — 1) 



/ 



where an empty entry means a zero entry. It has been established in many papers 
- see for example Mahmoud [17] , Chern and Hwang [6] or [1] and it can be easily 
checked that -Rg's (unitary) characteristic polynomial is 



XRgW 



n,, ,x . T(A + m) 



ml 



(2.3) 



where T denotes Euler's Gamma function. All eigenvalues are simple, 1 being the 
one having the largest real part. 



In the whole paper, X 2 will denote xr g ' s ro °t having the second largest real part 
and a positive imaginary part. 
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The famous phase transition on m-ary search trees already mentioned in the 
introduction is due to the fact that 

5R(A 2 ) > 1/2 if and only if m > 27. 

See for example pQ. The assumption 9ft(A 2 ) > 1/2 will be frequently used in the 
sequel. 

We adopt the following notations: 



VnGZ 



>0; 



z\ r(z + l) z{z-l)...{z-n + l) 



y n J n\T(z — n + 1) n\ 



2 + fc' 

Kk<m-1 



Ui(zi, . . . ,x m _i) = £ /cx fc ; 

l<fc<m-l 

U 2 («!,..., 2C m _i) = £ ( ^-1 ) 



l<fc<m-l 

A 2 + k - 1\ (2.4) 



l<fc<m-l 

1 / 1 



1>1 



Kfc<m-1 



1 / 1 

^2 



H m (X 2 ) \k(^ k )J i 



<k<m-l 



The linear forms u\ and u 2 are eigenvectors of G, namely G(u\) = u\ and G(u 2 ) = 
A 2 it 2 . The vectors t>i and v 2 are left eigenvectors of Rq, respectively associated 
with the eigenvalues 1 and A 2 . They satisfy Ui(vi) = u 2 {v 2 ) = 1 and U\{v 2 ) = 
^(vi) = 0. These eigendata had already been essentially computed in [I] and 
1201 . Note that, since A 2 is not real, u 2 and v 2 have nonreal coordinates. 



3 m-ary search trees and embedding 

In this section we present the connection between m-ary search trees and multi- 
type branching processes. This example of embedding of a discrete time process 
into a continuous time process has already been evoked by Janson in |12j . 

3.1 m-ary search trees 

We define here a discrete time Markov chain X DT = (X® T , n > 0) with values in 
N m_1 \ {0}. The z-th coordinate of X® T is denoted by Xn and has a "physical" 
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meaning detailed hereafter. The Markov chain X DT is a random walk defined by 
an initial vector X® T in N m_1 \ {0} and by the following transition probabilities: 
Vv G N™- 1 \ {0}, Vfc = 1, . . . , m - 1, 



kl k {v) 



q(v,v + w k ) = „ m _;\/ / - , (3.1; 



where the increment vectors w k are given in Section [2.11 and h(v) denotes the 
k-th coordinate of the vector v. 

Classically (see Norris [H] and for a synthetic exposition Bertoin [5]), this 
discrete time Markov chain is embedded in continuous time using a "Poissoniza- 
tion" of the time: given X DT , one can recover X = (X(t),t > 0) as follows. At 
time 0, X(0) = X DT . For any vector v G R m_1 , define^ 

m— 1 
k=l 

Let T\ be a random time exponentially distributed with parameter q(X® T ). For 
any time t G [0, T\[, let X(t) = X(0) = X^ T . At time n, X jumps from v = X(0) 
to v + Wfc with probability given by formula (13. 1 p . More generally, let r = and 
for any n > 1, define the n-th jumping time r n by 

n-l 

where are independent random variables having the same exponential distri- 
bution with parameter 1. Let 

X(t) = X(r n ) = X° T , Vt G [r n ,r n+1 [. 

At time T n +i, X jumps from v = X(r n ) to v + Wk with probability given by 
formula (13. ip . It is easy to see that this embedded process X(t) is the same one 
as the branching process defined in Section 12.11 

When X® T = (1, 0, . . . , 0), each Xn\i = 1, . . . , m — 1, can be seen as the 
number of nodes of type i in a tree T n : the sequence (T n , n > 0) is a sequence of 
random m-ary trees which grow by successive insertions of keys in their leaves. 
Each node of these trees contains at most m — 1 keys. Keys are i.i.d. random 



Note that q = u\ where u\ was defined by (|2.4 
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variables x^i > 1, with any diffusive distribution on the interval [0, 1]. The tree 
T n ,n > 0, is recursively defined as follows: T is reduced to an empty node-root; 
T\ is reduced to a node-root which contains Xi, T 2 is reduced to a node-root 
which contains X\ and X2, ■■■ , T m _i has a node-root containing X\ , . . . x m -\ . 
As soon as the (m — l)-th key is inserted in the root, m empty subtrees of the 
root are created, corresponding from left to right to the m ordered intervals 

h =]0,X(i) [,..., J TO =]Vi)' 1 t' where < x (i) < ••■ < £(m-i) < 1 are the 
ordered m — 1 first keys. Each following key x m , ... is recursively inserted in the 
subtree corresponding to the unique interval Ij to which it belongs. As soon as a 
node is saturated, m empty subtrees of this node are created. 

For each i = {1, . . . ,m — 1} and n > 1, X$ is the number of nodes in T n 
which contain % — 1 keys (and % gaps or free places) after insertion of the n-th 
key; such nodes are named nodes of type i. We don't worry about the number of 
saturated nodes. The vector X® T is called the composition vector of the m-ary 
search tree. It provides a model for the space requirement of the algorithm. One 
can refer to Mahmoud's book [17] for further details on search trees. 

Notice that, in this dynamics, the insertion of a new key is uniform on the 
gaps; it can be read on the transition probabilities (13. ip . 

3.2 Embedding 

The embedding properties are summarized in the following lemma. 
Lemma 3.4 

1) For any n > 1, the distribution of r n — r n _i is Sxp(n — 1 + Nq), where Nq 
is the number of free places in X(0): Nq = ui(X(0)). 

2) the processes (r n ) n >i and (A(r n )) n >i are independent. 

3) the processes (X(r n )) n >i and (X^ T ) n >i have the same distribution. 

Proof. Part 1) is a consequence of the fact that the minimum of k independent 
£xp( redistributed random variables is £xp(/c)-distributed, and that the total 
number of free places at time r n equals n — 1 + Nq. 

Part 2) is the classical independence between the jump chain and the jump 
times in such Markov processes. The initial states and evolution rules of both 
Markov chains in discrete time and in continuous time are the same ones, so that 
Part 3) holds. □ 

Convention . From now on, thanks to Part 3) of Lemma I3.4[ we will as usual 
suppose that the discrete-time process and the continuous-time process are built 



on the same probability space on which 



{X(r n )) n >i = {X? r ) n >i a.s.. (3.2) 

Remark . The important benefit we get with the embedding is the independence 
in the continuous-time process. This independence is the key point for the dislo- 
cation equations later on. 

4 Asymptotics and martingale connection 

In this section we present an order 2 expansion of the continuous time multitype 
branching process (X(t)) and its connection with the discrete time process (X® T ) 
defined in Section 13.11 

4.1 Asymptotics of the continuous time branching process 

With the notations of Section [2] and especially the formulae (12. 4p . the random 
vector X(t) admits the following order 2 expansion as t goes to infinity. 

Theorem 4.5 (Asymptotics of continuous time process) 

Suppose that m > 27. Then, as t tends to infinity, 

X(t)=e t £v 1 (l + e 1 (t))+2M(e x * t Wv 2 ) (l + e 2 (t)) + e 3 (t), (4.1) 

where 

• £ is a positive Gamma- distributed random variable with expectation N = 
ui(X(0)) (total weighted number of particles at time 0), 

• W is a complex-valued random variable that admits moments of all orders p > 1 
and whose expectation equals U2(X(0)), 

• the real-valued random variables Ei(t) and £2{t) tend to as t tends to +oo, 
almost surely and in any L p -space, p > 1, 

• the random vector £s(t) is o (e A2i ) as t tends to +oo, almost surely and in any 
L p -space, p>l. 



In the whole paper, W denotes our hero, namely the limit complex-valued ran- 
dom variable of the second order term in X(t) 's expansion, as in Theorem \4-5[ 
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Remark 4.6 One can reformulate fl^.lp as follows: 
in any basis of the form (v\, 9ft(t> 2 ), S(f 2 ), . . . ), 

• X(t) 's first coordinate has the expansion e*£ + o(e l ), 

• X(t)'s component on Span R (^(i^), ^(^2)) has the expansion 23? (e X2t Wv 2 ) + 
o(e X2t ), 

• all other coordinates are o(e x ' 2t ). 

Proof of Theorem H3J Denote A the endomorphism of IR" 1 " 1 having l R G as 
matrix in the canonical basis. Let also M(t) = exp(— tA)X(t), for any t > 0. 
By standard arguments from multitype branching process theory, (M(t))t>o is a 
vector- valued martingale. Since m > 27, the real part of A2 belongs to ]l/2, 1[ 
so that the projected martingales u\(M(t)) and u 2 (M(t)) converge in L p for any 
p > 1. For proofs of these results, see for example Athreya and Ney [T] or Janson 
[T2] (especially Lemma 10.2 of Janson's paper for the L p -boundedness, X being 
here an irreducible process in the sense of p2]). The random variables £ and W 
are respectively defined by 

£ = lim ui (e- tA X(t)) = lim e~V (X(t)) , 

i— >+oo t— >+oo 

(4.2) 

W = lim u 2 (e~ tA X(t)) = lim (X(t)) . 

An alternative proof of the L p convergence can be made using the techniques of 
[2Tj , as developed in [5] for two-colour urn processes. In particular, £'s distribution 
is attained by explicit computation of its moments: for any nonnegative integer p, 
an elementary computation shows directly from (12.21) that the (so-called reduced) 
polynomial 

Q := Mi (iti + 1) (1*1 + 2) . . . (ui + p - 1) 

is an eigenvector for X's infinitesimal generator G, associated with the eigen- 
value p. Thus EQ(X(t)) = e pt Q(X(0)) for any t. Besides, because of 
Q(X(t)) = e pt ^ p {l + o(l)) as t tends to infinity, almost surely and in L 1 . Finally, 
the last two equalities provide 

W = Q(X(0)) = r( ^ 0+p) . 
^ ^ v ^ n r (No) 

This shows that the law of £ is a Gamma distribution with parameter Nq since a 
Gamma distribution is completely determined by its moments. The matrix Rq 
is diagonalizable on C since all roots of its characteristic polynomial are simple 
(see (12. 3p ). Extending notations (I2.4p . let (ma)agSp(^) be a basis of linear forms, 
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each u\ being an eigenform of G associated with the (complex) eigenvalue A. Let 
also (i^AeSpiU) be the dual basis of («a)agSp(^i), eac h v x being thus a vector that 
satisfies u x (v^ L ) = 5 X)fJ- (Kronecker's notation). Note that one can choose u\ 2 = u 2 
and, consequently, v\ 2 = v 2 (cf. notations (12.4)) ) . 

For any t > 0, split the spectral decomposition of the vector X(t) with respect 
to G into four terms: 

x(t) = MX(t)).v x = X^t) + X 2 (t) + X 3 (t) + X 4 (t), 

where 

' X 1 (t) = u 1 (X(t))v 1 , 
X 2 (t) = u X2 (X(t))v X2 + u- 2 {X{t))i % , 

x 3 (t)= «A(*(t)K 

l/2<5ftA<5RA 2 
k = J2nX<l/2M X (t))vX- 

Note that this partition of Sp(.A) is valid because \ is not an eigenvalue of A 
as can be checked from (j2.3p . We deal separately with these four components of 
X(t). Define s 3 by e 3 (i) = X 3 {t) + X 4 (t), for any t > 0. 

• The formulae ( 14. 2 j) provide directly the asymptotics 
f X 1 (t) = (e*e + o(e*)K, 



[ X 2 (t) = (e A2f W + o(e A2 *))w 2 + (e«f + o(e x ^))v 2 = 29ft ((e A2 'I¥ + o{e x ^))v 2 ) , 

leading to the first two terms of the expansion (14. ip . 

• Suppose that A is an eigenvalue of A such that ~ < 3ft A < 9ftA 2 . Then, with 
the same general arguments as in the very beginning of the proof, it can be seen 
that 

u x {M{t)) = e~ tx u x (X(t)) 

and that (u\(M(t)))t>o is a convergent martingale, bounded in any L p , p > 1. In 
particular, w^(X(t)) = o(e A2 *) as t tends to infinity, almost surely and in any L p , 
p > 1. This shows that ^(t) is o(e A2 *) when t — >■ +oo. 

• It remains to deal with the small eigenvalues, namely with all A such that 
3ftA < \. 
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Lemma 4.7 Suppose that A is an eigenvalue such that 3ft A < ~ and let r\ > 0. 
Then, e~^ +r, ^u x (X (t)) zs bounded almost surely and in any L p -space, p > 1. 
The proof of this lemma is given just hereafter. Therefore, if 9ftA < |, then 

~ A2 'u A (X(t)) = e (1/2+ "- A2)t \e-^ +ri)t u x (X(t)) 



e 



— y 

t— >oo 



almost surely as soon as < ^ < 3ftA 2 — |. Such r\ exist because 3ftA 2 > \. This 
shows that X 4 (t) is o(e A2< ) when t — > +oo. The same argument holds for the L p 
convergence, making the proof complete. □ 

Proof of Lemma 14.71 The main idea consists in taking advantage of the 
following fact: when t belongs to the interval [r n ,r n +i[, the vector X(t) remains 
equal to X® T . This being considered, we make use of the moment bounds of the 
discrete time process that can be found in [21] (Theorem 3.4 (1)): when 3ftA < |, 

Vp>l,Ve>0, E\u x {X° T )\ p = O (n p ^ +£) Y n -)■ +oo. (4.3) 
• Almost sure bound: we prove that 

lim P (3t > 0, e~^ +Tl)t \u x (X(t))\ > C] = 0, (4.4) 

which suffices to get the almost sure boundedness. Let C > 0, rj > and let A be 
an eigenvalue such that 3ftA < \. The jump time r n tends almost surely to +oo 
which is a classical result that can be deduced from Lemma I3.4[ so that 

p (at > 0, e-^MX(t))| > 6) < { 3t e K,r n+1 [, e-^\u x (X(t))\ > c) 

n>0 

Since X(t) = X® T for any t G [r n , r n+ i[, this leads to 

P (at > 0, e~^ +r ^\u x (X(t))\ > c) < J2 ¥ (^,f T )| > Ce^+^A 

n>0 

Conditioning with respect to r n , using Markov inequality and the fact that r n 
and X® T are independent, one gets successively, for any p > 1: 



P (3t > 0, e-d+^|u A (X(t))| > c) < (V (Ju A (X r ? T )| > Ce ( ^ 



' r 



n>0 
n>0 



Eh (A^T 



(7p e p(|+'?) T ™ 



n>0 
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The density of the n-th jump time r n is the function 

n£l i — > ne~ u (1 - e"")™" 1 l R+ (u), 
so that its Laplace transform can be elementarily computed: for any s > 0, 

n\T(s + l) 



E(e 



r(s + l)n s , n ->■ +00. 



r(s + l + n) 

Together with this leads to: V77 > 0,V£: > 0, Vp > 1, 

EK(inr«H*)=°(^)' 



rt — > +00 



which is the general term of a convergent series as soon as one takes e < 77 and 
p > -3-. Finally, letting C tend to infinity shows (14. 4p . 

• Bound in L p -space: let p > 1 and t > 0. Then, 

-^ +r?) 'u A (X(t)) " = e - ( 5+"^E |u A (X(t))| p . 
p 

Using the relation with the discrete time process (X^ T ) n , one has successively 



<h+^ Ux (x(t)) 



e y 2 



n>0 



l T n <t<T n + 



i|«A (*(*)) 



n>0 



e -(5+?7)pt 



^E(l Tn < t<rn+1 )E(|n A (X n DT )| p 



n>0 



where the last equality holds due to the independence between r n and X^ T . 
Besides, r n and r n+ i — r n are independent and r n+ i— r n is + A r )-distributed 

(see (I3.4p ). so that, using the density of r n written above, one gets 



^(lr n <t<T, l+ i 



E flt> T „E ^l r „ + i-r„>t-r„ I 7"n 

E(l t > rn e-( n+7V °)('- r ")) 







< ne 



-(n+l)t 



f\e u -l) n - l e u du= (1-e-*) 



n n — t 
e . 
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Thus, 



P < e -t e -{\+i)pt- 



j2(i-e~r®{\^(Xn T )\ 



n>0 



Let now e > 0. On one hand, (14.31) implies that 

E(|«A (Xf r )| p )=o(n^). 



On the other hand, Stirling's formula applied to generalized binomial coefficients 
yields classically that for any a G C, 



[z»](l-z) 



-a-l 



n a ( ( 1 

r(^( 1 + U 



where the notation [z n ]A(2i) means the coefficient of z n in the power expansion 
of A(z) at the origin. Consequently, 

e(K W T )r) = o(H(i-,r 1 -^)). 

This implies that for any 5 > 0, there exists a constant C e such that for any t > 0, 

-(^V (X(t)) " < C B eV3^ (1 - (1 - e - t ))~ 1 " p( ^ +£) = C £ e- pt ^- £ \ 
p 

It suffices to take e = r]/2 to conclude that the L p -norm of e~ ( -2 +r i) t ux (X(t)) is 
bounded above. □ 

Remark 4.8 The distribution of W is infinitely divisible, because it is the limit 
of infinitely divisible ones, obtained by scaling and projection of infinitely divisible 
ones. Indeed, in finite time, for any xq G M m_1 , denote by (X Xo (t),t > 0) the 
process (X(t),t > 0) defined in Section \2J\ starting from initial state xq. By the 
branching property 

X X0 (t) £ [n]Xa(t), 

where the notation [n]X denotes the sum of n independent copies of the random 
variable X. This fact has already been noticed by Janson (fTEj . proof of Theo- 
rem 3.9). 
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4.2 Martingale connection 



In this subsection, we use the embedding equality (13. 2 p to deduce connections 
between the asymptotic behaviours of X® T when n — > +oo and X(t) when t — > 
+oo. 

For m > 27, it has been proved in [I] and [20] that 

X° T = nv x + $t(n X2 W DT v 2 ) + o(n U2 ) a.s. and in L p , Vp > 1, (4.5) 

where v± , v 2 are the deterministic vectors defined in (12 .4p , W DT is a complex- 
valued martingale limit, and the moments of W DT can be recursively calculated. 

Proposition 4.9 The following two assertions hold: 

lim ne~ Tn = £ a.s. and in L p , Vp > 1, (4-6) 

n— >+oo 

W = £ A2 W DT a.s. with £ and W DT independent. (4.7) 

The equality ( 14. 7p is commonly referred to as "martingale connection" . 

PROOF. We first prove ( 14. 6p . Applying the first projection to the embedding 
equality ( 13. 2p . we obtain that 

uipT(O) = Ul (X° T )) a.s., 

where U\ has been defined in (12 ,4p . This is the total number of free places at time 
r n , and is equal to n — 1 + N = n(l + o(l)). Therefore, by ( I4.2p and the fact that 
the splitting times r n tend almost surely to +oo when n goes to +oo, we have 

^ = lim e _i Mi(X(t)) = lim ne~ Tn a.s.. 

t— >+oo n— >+oo 

This gives ( 14. 6p . 

We then prove (14. 7p . Applying the second projection to the embedding equal- 
ity (13. 2p we obtain 

u 2 (X(r n )) = u 2 {X° T )) a.s., 

where u 2 has been defined in ( 12 .4p . Using again (14. 2p and the fact that r n goes 
to +oo when n goes to +oo, we get 

W = lim e- X2t u 2 (X(t)) = lim e- X2T "u 2 (X° T ). 

Therefore (14. 7p follows from (14. 6 p and the asymptotics in discrete time given in 
( 143j) . □ 
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Remark 4.10 Fill and Kapur (f^j) proved that W DT is the unique solution in 
the space of probability distributions with a given mean and finite second absolute 
moment of the fixed point equation 

m 

Z = ^(V fe ) A2 ^ W , (4.8) 

k=l 

where the Vk are the spaces in the statistical order of (m — 1) i.i.d. random 
variables uniformly distributed on [0,1]. Because M.-valued stable distributions 
are solutions of the fixed point equation ^4-8j) when X is a real number, it is 
somehow natural to ask whether a X-stable distribution is a solution of Eq. H^W 
for a complex number X. By X-stable we mean operator-stable when the operator 

is given by a two dimensional matrix X = a + ir = ^ ^ T ^ as introduced by 

Sharpe J^] /. It is known since Hudson et al. fI7]/ that a X-stable distribution has 
infinite moments of order p for p > 1/3?(A). Consequently, neither W nor W DT 
(which have moments of any order) can be stable distributions. 



5 A distributional equation 

In this section we derive a distributional equation satisfied by the limit variable 
of the continuous-time branching process with an appropriate norming. We shall 
see that this equation characterizes the limit distribution. 



5.1 Vectorial finite time dislocation equations 

Let us write dislocation equations for the continuous-time branching process at 
finite time t. We write Xj(t) for X(t) when the process starts from X(0) = ej, 
where e,j denotes the j-th vector of the canonical basis of IR m_1 (whose j-th 
component is 1 and all the others are 0). This means that the process starts from 
an ancestor of type j. 

Notice that the first splitting time t\ changes of distribution depending on 
the ancestor's type; denote by T^\,j = 1, . . . , m — 1, the first splitting time when 
the process starts from X(0) = e,-. Thus ry) is Exp{j) distributed. 

The branching property applied at the first splitting time gives: 
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Vi > Tu < 



Xx(t) =X 2 (t-r {1) ), 

X 2 (t) ±X 3 {t-T {2) ), 



X m -2(t) — X m _i(t — T( m _2), 



(5.1) 



X m _i(t) = [m]Xi(t - r (m _i)), 

where the notation [m]X denotes the sum of m independent copies of the random 
variable X. 

In the following, we denote 

T = t {1) ^ hT (m _i), (5.2) 

where the tu) are independent of each other and each ry) is Sxp(j) distributed. 
Let us give some elementary properties of T that we shall need. It is easy to see 
that T has density 



f T (u) = (m - l)e-«(l - e"«) 



u\m—2 



1r+(m), «e 



(5.3) 



so that e~ T has a Beta distribution with parameters 1 and m — 1. A straightfor- 
ward change of variable under the integral shows that for any complex number A 
such that 3R(A) > 0, 



Ee 



-AT 



r+ao 

I e- Xu f T {u)du = (m-l)B{l + \m-l) 
Jo 

(m- 1)! 



nr=i 1 (A+^)' 

where i? denotes Euler's Beta function: 



(5.4) 
(5.5) 



B(x,y) 
In particular, 



mE\e 



-XT I 



iw-irfu = r ( g ) r ^) sft x > o, % > 0. 
T(x + y) 

< 1 if ft(A) > 1, 
= 1 if 3?(A) = 1, 
> 1 if ft(A) < 1. 



(5.6) 
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5.2 Distributional equation satisfied by the limit variable 

After projections of variables Xj(t) (the process starts from X(0) = ej) with 
scaling with e~ X2t and taking the limit when t goes to infinity, we get the variables 



W. 



j ■' 



lim e- A2 *M 2 (X j (t)), 



t—t+oo 



so that the system 05. in on Xj(t) leads to the following system of distributional 
equations on Wf 



2- 



W 2 = e~ X2T ^W 3 , 



< ••• 



(5.7) 



W, 



m-2 



w, 



m—1 



- A 2 T (m-i)[ m ]Wi. 



This shows that W\ is a solution of the following fixed point equation: 

Z k e" A2T (Z (1) + --- + Z (m) ), (5.J 

where are independent copies of Z, which are also independent of T. 
In terms of the Fourier transform 

(p{t) := Eexp{i(t,Z}} = E exp{«3?(tZ)}, t G C, 

where (x, y) = = 9ft(x)9ft(y) + S(x)S(y), the equation (15.81) reads 



ip m (te- X2U )f T (u)du, t E C, 



(5.9) 



where /y is defined by (15. 3p . Notice that this functional equation can also be 
written in a convolution form: if $(t) := ip(e X2t ) for any t G C, then $ satisfies 
the following functional equation: 



$(*) 



$ m (t - u)f T {u)du, t e C. 



(5.10) 



In the following sections, we prove that the distributional equation (15.81) char- 
acterizes the law of W\ and we get several results on W\\ for example W\ has 
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a density on the whole complex plane, and admits exponential moments. All 
these results appear as a particular case of a slightly more general situation given 
hereafter. From now on, for any complex number A, consider the distributional 
equation 

Z = e" AT (Z (1) + --- + zW), (5.11) 

where Z^> are independent copies of Z, which are also independent of T. On 
Fourier transforms, it reads 

i>+co _ 

(p(t) = / <p m (te- Xu )f T (u)d Uj t G C, (5.12) 
Jo 

where fr is defined by ( 15. 3j) . 

Notice that when Z is a solution of the distributional equation (15. lip , with 
finite and non zero first moment, then A is a root of the polynomial function (12. 3p . 
In particular, A is an algebraic number. 



6 The smoothing transformation 

A solution of the distributional equation (I5.1ip is a fixed point of the associated 
smoothing transformation defined hereafter by (16.11) . Endowing a suitable space 
of probability measures with two distances, we prove that the smoothing transfor- 
mation is a contraction for both metrics. This provides two alternative approaches 
for the study of Eq. (15. lip by contraction method. Using the Wasserstein distance 
as a first metric, we adapt the classical contraction method developed in [10] , [22] 
and [23]. The second metric is defined in terms of Fourier transforms of measures; 
it provides a short proof of our result. 

For any complex number C, let M. 2 {C) be the space of probability distribu- 
tions on C admitting a second absolute moment and having C as expectation. 
Let A be a complex number. For any probability measure \x on C, let 

Kfi := £(e" AT (ZW + • • • + Z< m) )), (6.1) 

where T is given by (15. 2p . are independent random variables of law fj,, which 
are also independent of T. Following Durrett and Liggett [8] who considered the 
case of real random variables, we call K the smoothing transformation. Note that 
K depends on m and A. 

Lemma 6.11 If X is a root of the characteristic polynomial 112. 3\) such that 
3?(A) > — \ and if C is any complex number, then K maps A^C) into itself. 



Proof. Since dt(X) > —1, the random variable e AT has an expectation. Further- 
more, by (JE3D, mEe~ AT = 1 as A is a root of (J23]). This ensures the conservation 
of the expectation by K. Since 5J(A) > — |, then E|e~ AT | 2 < oo and K\i admits 
a second absolute moment whenever \i does. Therefore Kfi G M. 2 {C) whenever 
fieM 2 (C). □ 

Notice that a solution of Eq. (15. lip is a fixed point of K. We shall use the 
Banach fixed point theorem for two different metrics on Ai 2 (C) to study the 
existence and uniqueness of solutions of Eq. (15. lip . 



6.1 Wassertein distance 

Let d 2 be the Wasserstein distance on M. 2 {C) (see for instance Dudley [7]): for 
//,*/ G M 2 (C), 

d 2 (a,u) = ( minE(|X-F| 2 ) 

\(X,Y) V ' 

where the minimum is taken over couples of random variables (X, Y) having re- 
spective marginal distributions fi and v\ the minimum is attained by the Kantorovich- 
Rubinstein Theorem - see for instance Dudley [TJ, p. 421. With this distance d 2 , 
A4 2 (C) is a complete metric space. 

Theorem 6.12 Let A G C be a root of the characteristic polynomial A2.3\) such 
that SR(A) > \, and let C G C. Then K is a contraction on the complete metric 
space {M. 2 {C),d 2 ), and the fixed point equation A5.11\) has a unique solution Z 
mM 2 (C). 

We now come back to the limit variable W± of m-ary search trees. Since 
EWi = 1 and IE | PVi | 2 < oo, the following corollary is a direct consequence of 
Theorem 16.121 applied for A = A 2 . 

Corollary 6.13 The distribution of the limit complex random variable W\ is the 
unique solution in the space J\A. 2 {1) of the fixed point equation Ii5.11\) . 

Proof of Theorem 16. 121 We argue as in [10] , [22] and [23] where real random 
variables were considered. 

By the Banach fixed point theorem, it suffices to show the contraction prop- 
erty. Let fi, v G M. 2 {C). Let (X,Y) be a couple of complex- valued random 
variables such that C(X) = ji, C(Y) = v and d 2 {\i,v) = a/E|X — Y\ 2 . Let 
(Xi, Yi), i — 1, . . . , m be m independent copies of the ^-optimal couple (X, Y), 
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and T be a real random variable with density fa defined by ( 15. 3p . independent 
from all (Xi,Yi). Then, 

m m 

C{e~ XT X i) = K V and C ( e ~ XT Y ^> = Kv > 

i=l i=l 

so that 

d 2 (K/i : Kvf < E '~~ A: 



i=i 



i=i 



E 



-AT 



i=l 



E |e~ AT rE 



E e 



-ATI 



,i=l 



= mE |e | d 2 (A*, ^) • 

Since 23?(A) > 1, we have mE |e _2AT | < 1 (see (15.61) ). Therefore K is a contraction 
on A^2(C) and the proof is complete. 

□ 



6.2 Distance defined with Fourier transforms 

We now give an alternative approach for the characterization of the limit dis- 
tribution via Fourier analysis. We define another distance d\ on M. 2 {C) as 
follows. Take fi, v £ A4 2 (C) and denote respectively tp and ip their character- 
istic functions. By definition of Ai 2 (C), both p> and ip admit the expansion 
(p(t) = 1 + i{t, C) + O (\t\ 2 ) when t tends to 0. Therefore, one can define d 2 (fi, v) 
by 

x \<p(t)-m\ 

d 2 ((i,u)= sup -— . 

t€C\{0} \t\ 

Clearly, d 2 ((p,ip) < oo, and d 2 is a distance on Ai 2 (C). It can be easily checked 
that (Ai 2 (C),d 2 ) is a complete metric space. 
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The following result is a counterpart of Theorem 16.121 It gives an alternative 
proof for the existence and uniqueness of the solution of Eq. (15.111) in the class 
of probability measures on C with a given mean and finite second moments. 

Theorem 6.14 Let A G C be a root of the characteristic polynomial Ii2.3\) such 
that 9fc(A) > \, and let C G C. Then K is a contraction on the complete metric 
space (A^2(C) , d 2 ) , and the fixed point equation Ii5.11\) has a unique solution Z 
inM 2 (C). 

Proof. Thanks to Banach fixed point theorem, it suffices to prove that K is a 
contraction on M. 2 {C) equipped with the metric d 2 . Let fi, v G M. 2 {C) and let 
if and ip be their respective characteristic functions. An elementary computation 
shows that the Fourier transform of K[i is t \— > Eip m (e~ XT t) with a corresponding 
formula for v. We have |</?| < 1, \i\>\ < 1, so that 

E\ip m (e- jT t) -i; m (e- jT t)\ < mE\ip(e- jT t) -?P(e- jT t)\. 

Together with the inequality \<f(z) — ip(z)\ < d^ifi, v)\z\ 2 applied to z = e" XT t, 
this implies that 

d* 2 (K/i,Ku) < mE(e- 2diWT )d* 2 {fi,u). 

Since 23?(A) > 1, we have mE ^e _25R( - A ^ T ) < 1 (see (I5.6P ). Therefore the above 
inequality shows that K is a contraction on (A^C); d 2 ). □ 

Remark 6.15 Denote ^(C) the space of Fourier transforms of elements of 
A4 2 (C). When X is a root of the characteristic polynomial $2.3\) such that 3?(A) > 
2, the smoothing transformation K can be identified as a map (also denoted by K) 
on .7-2 (C) given by 

(K<p)(t) :=Ep m (e- jT t), t G C. (6.2) 

The proof of Theorem \6.14\ a ^ so shows that K is a contraction of .7-2 (C) for the 
metric (also denoted by d 2 ) defined on ^(C) by 

d> 2 M):=suv Mt \- m . (6.3) 

Let T> 2 (C) be the space of all continuous functions (f : C — > C that admit an 
expansion ip{t) = 1 +i(t,C) + 0(|t| 2 ) at and such that ||y?||oo < 1- Clearly, 
T> 2 (C) contains J- 2 (C). One can show that formula Ii6.2\) defines a mapping from 
T> 2 (C) into itself and that K is a contraction for the metrics defined by (TjOj). 
This provides a proof of existence and unicity of solutions of Ii5.12\) on T> 2 (C). 
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Remark 6.16 One can deal with the convolution equation A5.10\) by arguments 
of the same vein. Similar computations show that this equation has a unique 
solution in the space 8%{\ C) of continuous functions $ : C — > C that admit an 
expansion $(z) = 1 + i(e Xz ,C) + O (e 2Az ) when \z\ tends to +00 and such that 
Halloo < 1- This result appears one again as a consequence of Banach Theorem 
on £2 (A, C) for the metric 



d($, = sup 

zee 



e 2Xz 



As a consequence, this shows in particular that the Fourier (complex) transform 
<p of W\ satisfies: for any wgC and for any determination of the logarithm, 



A 



2 



where $ is the unique solution in £2(^2, 1) of Eq. A5.10\) . This result is the 
reversed version of the change of variable §(z) = ip(e X2Z ) that led from A5.9\) 
to tFM) . 

7 Density and support 

In this section we prove results on the absolute continuity and on the support 
of solutions of the distributional equation (15. lip via Fourier analysis. As appli- 
cations, we show that the distribution of the limit variable W\ of the multitype 
branching process (X(t)) has always a density and that its support is the whole 
complex plane. 

Theorem 7.17 Let X be a complex number such that A ^ 1 and 9ft(A) > 0. Let Z 
be a complex-valued random variable solution of the distributional equation Ii5.11\) 

Z ± e - XT {zW + --- + zW), 
with K\Z\ < 00 and KZ ^ 0. Then the following assertions hold: 
(i) the support of Z is the whole complex plane C; 

(11) as \t\ 00, Ee i{t ' z) = 0(\t\- a ), for any a e]0, ^7tt[; 

(Hi) the distribution of Z has a density with respect to the Lebesgue measure 
on C. 
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Remark 7.18 When X = 1, the distributional equation A5.11\) becomes 



X = e- T (X« + ... + XM). (7.1) 

By Section® it admits a unique solution in the space A^C) of probability mea- 
sures on C, with a given mean C . Moreover, from the dislocation equations 115. 
a similar argument shows that 

£:= lim e~WX(t)) 

t— >+oo 



is a solution of this equation. By Theorem 4-5, £ is Gamma- distributed. Therefore 
the unique solution of J7.jp in M. 2 {C) is C7 where 7 is Gamma(l)- distributed, 
and its support is the half line CR + . 

The following corollary gives the main result for the limit variable W\ of the 
multitype branching process. It is a direct consequence of Theorem 17.171 since 

mv l = 1. 

Corollary 7.19 The distribution of W\ admits a density with respect to the 
Lebesgue measure on C, and its support is the whole complex plane C. More- 
over, as \t\ -»■ 00, Ee i{t ' Wl) = 0{\t\- a ) for each a e]0, [■ 

The proof of Theorem 17. 171 runs along the following lines. Let ip be the Fourier 
transform of any solution Z of ( 15. IIP . We prove that tp is in L 2 (C) because it is 
dominated by \t\~ s for some 5 > 1 so that the inverse Fourier-Plancherel transform 
provides a square integrable density for Z. The guiding idea consists in adapting 
usual methods (developed in [H] and [15]) used for positive real- valued random 
variables to our complex-valued case thanks to the function defined for r > by 

ip(r) = vMx\<p(t)\. 

\t\=r 

From now on, 

A = e~ XT . 
We proceed by a series of lemmas. 

The first lemma concerns a property of the support of Z. For a complex- valued 
random variable Z and for a complex number z, by definition, 

z E Supp(Z) > 0, F(\Z - z\ < e) > 0. 
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Lemma 7.20 Let z G C. Then 



z G Supp(Z) D(0, \z\) C Supp(Z), 

where D(0, \z\) denotes the open disc with center and radius \z\. 

PROOF. We first prove the following implication: 

[a G Supp(A) and z G Supp(Z)] =>- maz G Supp(Z). 

Indeed, let e > 0, a G Supp(A) and z G Supp(Z). Let also . . . , 

be independent copies of Z. Then, with positive probability, \A — a\ < e and 

\Z^> — z\ < e for any k. Therefore, with positive probability, 



\A(Z<U + • ■ ■ + ZW] 



maz\ 



mz(A - a) + Aj2( z(k) ~ z ) 



< me\z\ + (|a| + e) me. 



The positive real e being arbitrary, this shows that maz G SuppA(Z« + •■■ + 
Z^) which implies that maz G Supp(Z) by Eq. (15. lip , proving the claim. 

Let z G Supp(Z). Since Supp(T) = M + (see (I5.3P ). the claim implies that for 
any t > 0, me~ xt z G Supp(Z). Iterating this property of Supp(Z) shows that 

{m n e- xt z, n G N, t G 1R+} C Supp(Z). (7.2) 

Since the support of a probability measure is a closed set, to show that D(0, \z\) C 
Supp(Z) it suffices to prove that [m n e~ xt , n G N, t G M + } is everywhere dense 
in the unit disc. Taking logarithm, we show hereunder that 

G := Nlogm + 2?7rN - XR + 

is everywhere dense in the half-strip 

B := {x + iy, x < 0, - 2ir < y < 0} 

which implies the desired result. 

Let a and r denote respectively the real and imaginary parts of A. Remember 
that, as soon as Z is a solution of Eq. (l5.11|) such that E|Z| < oo and EZ 7^ 0, 
then A is a root of xr g defined by (12. 3p ; this implies that A is an algebraic number, 
and so is a jr. 
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Let us prove that p := in/ logm is a transcendental number. In fact, if p were 
algebraic, then by the Gelfond-Schneider theorem^, m p would be a transcendental 
number; but this is impossible because m p = exp(iir) = —1. Therefore p is a 
transcendental number. 

It follows that 27rcr/(r logm) is not an algebraic number, hence not a rational 
number. So by a classical result, 

N log m — N is a dense subset of R. 

T 

Let b = x + iy G B with x < and — 2n < y < 0, and let e > 0. Approximating 
the real number x — — by an element of Nlogm — N^ 21 , we take n,k G N such 
that 



(n log m — k ) — (x — — ) 



< e. 



Therefore 



2&7T — y 

o — g\ < £ ) where (7 = nlogm + zzkn — t A G Q with t = > 0. 

T 

This completes the proof of Lemma 17.201 □ 

Lemma [7.201 leads to the following key property of ip, which will imply that 
the characteristic function (p of Z satisfies the Cramer condition and then the 
Riemann-Lebesgue condition. 

Lemma 7.21 Vr > 0,ip(r) < 1 . 

Proof. Obviously, ^(0) = 1 and ip(r) < 1 for any r > 0. Suppose that r > 
is such that ip(r ) = 1. Take thus z G C and 9 G E such that 



Pol 



r and Ee i{z °' z) = e 



i0 



The complex random variable e l (( z o, z )- e o) j s Q f m ean 1 and takes its values on the 
unit disc, so that it is almost surely equal to 1. This implies that Supp(Z) is 
contained in a set of countably many parallel lines of the complex plane. This 
contradicts Lemma 17.201 since such a set of lines is negligible with respect to the 
Lebesgue measure on C. □ 



3 The Gelfond-Schneider theorem states that if a and b are algebraic numbers with o ^ 0,1 
and if b is not a rational number, then any value of a b = exp(6 log a) is a transcendental number. 
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Remark 7.22 The preceding arguments show the following assertion: for any 
complex-valued random variable Z , if |E(e^ 2 °' z ^)| = 1 for some z G C\ {0}, then 
Supp(Z) C a + bZ + cR for some a,b,c G C (a set of countably many parallel 
lines). The algebraicity of A that leads to the proof of Lemma \7.20\ can thus be 
seen as a nonlattice assumption on the fixed point equation A5.11\) . 

We now prove that the characteristic function <p satisfies the Riemann-Lesbesgue 
condition. 

Lemma 7.23 lim ip(r) = 0. 

r— >+oo 

Proof. We argue as in the proof of Theorem 3.1 of [2] or Lemma 3.1 of |15j . 
Notice that from the distributional equation (15.121) we have 

V>(r) < E(^ m (r\A\)). (7.3) 

• We first prove that lim sup ip(r) = or 1. By Fatou's lemma, 

lim sup ip(r) < Elimsup , ?/' m (r|A|) = I lim sup ip(r) 

r— >+oo r— >+oo \ r— >+oo 

the last equality coming from P(|A| > 0) = 1. So the real number I := lim sup r ip(r) 
satisfies both I < 1 and I <l m \ this implies that / = or 1. 

• Suppose that limsup r ^(r) = 1. By Lemma XI if}(l) < 1. For any e G 
]0,1- ^(1)[, define 

{r\{e) = max{r g]0, 1[, ip{r) = 1 — e}, 
r 2 (e) = min{r > 1, ip(r) = 1 — e}. 

These quantities are well defined because ^(0) = 1 and i/j is continuous. Then 
^(ri(e)) = ip{r2{e)) = 1 — e and for any r G [ri(e), ^(e)], i/>(r) < 1 — e. 

Let us prove that r\{e) goes to when e tends to 0. Take any limit point p 
of Txie). Since ip is continuous, ip(p) = 1 — e which implies by Lemma [7.211 that 
p = 0: the only possible limit point is 0. 

By (17.31) . we have 

^(r) < E^(r|A|). 
Iterating this inequality we see that for all n > 1, 
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where (|A|)i>i are independent copies of \A\. With the notation 

X n (r,e) :=¥(r 1 {e)<r\A 1 \...\A n \ <r 2 (e)), 
we have for any r > 

ip(r) < (1 - e)A n (r, e) + 1 - A n (r, e) = 1 - eA n (r, e). 
Again by (15321) . 

1 _ e = ^,( r2 ( e )) < E0 m (r 2 (e)|A|) < E (l - £ A n (r 2 ( £ )|A|, e) 
In other words 



E(l- (l-£A n (r 2 (£)LA|, £ 



< 1. (7.4) 
We are going to pass to the limit in the above ratio when e tends to 0. Rewrite 

X n (r 2 (e)\A\,e) = P < \Al\Ai\ . . . \A n \ < , 

and remember that r\{e) < 1 < r 2 (e) and that r\{e) goes to when e tends to 0, 
so that < nM 

goes to when e tends to 0. Consequently, 

X n (r 2 (e)\A\,e S j — j p(o < . . . \A n \ < l) = a.s., 

where, for any x > 0, 

:=P^|A|...|A n | < l). 



Therefore 

1- (l-e\ n (r 2 (e)\A\,e 



m/i n (|A|) a.s.. 



e e^O 

The above ratio is a function of e, uniformly bounded on the compact set [0, 1 
"0(1)]) so that by dominated convergence and (17. 4p . 



E(l- (l-e\ n (r 2 {e)\Ale 



Besides, by Markov inequality 

1 - /in 0*0 < xEQAxl . . . \A n \) = x (E\A\) 



mEfi n (\A\) < 1. (7.5) 
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Since 3ft(A) > 0, E\A\ =E\e~ XT \ < 1 (see (J53D), which implies that lim^-^ fi n (x) = 
1, so that 

lim E/i n (\A\) = 1 

n— >oo 

by dominated convergence. This contradicts (17. 5p because m > 2. □ 

We shall need an information about the decay rate of <p(t), of the form tp(t) = 
0(\t\~ 5 ) for some 5 > when |t| — > oo. To this end, we shall use the following 
Gronwall-type technical Lemma of [H] (see also Lemma 3.2 in |15j). 

Lemma 7.24 [Lffl Let ip : R + — > R + be a bounded function and let B be a positive 
random variable such that for some constants p g]0, 1[, a > 0, C > and for all 
r > 0, 

ip{r) < pEip(Br) + Cr~ a . 
If pE(B~ a ) < 1, then ip(r) = 0(r~ a ) as r -> oo. 

This is Lemma 4.1 of [H]. It can be proved as follows. Let {B{\ be independent 
copies of B. Then by induction, for all n > 1 and all r > 0, 

^(r) < p n Eip(B 1 ...B n r) + CV~ a [l + pE(£?~ a ) + ... + (pE(B- a )) n - 1 ]. 

Letting n — > oo we see that for all r > 0, 

ip(r) < Cr' a /[1 -pE(B- a )}. 

Lemma 7.25 For all a g]0, [, as r ^ oo, 

ux( A) 

^(r) = 0(r~ a ). 

Proof. We have already seen from the distributional equation (I5.12p that 

^(r) < E(^ m (r\A\)), 

where A = e~ XT . By Lemma I7.23[ for any e > 0, there is some r £ > such that 
Vr > r e , ^(r) < e. So 

^(r) < e m_1 E^(r|A|) +P(r|A| < r e ). 



Therefore by Markov inequality, for a G]0, 




ip{r) < e m - l E4){r\A\) + r~ a (r £ ) a E(|A|" a ). 
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By (JE3D, E (l^l" a ) = (m- 1)5(1 - aU(X),m - 1) < oo. Taking e > small 
enough such that £ m-1 E(|A| _a ) < 1, we see that the desired result follows from 
Lemma [7.241 and the preceding inequality on ip(r). □ 

We can now finish the proof of Theorem 17.171 

Proof of Theorem 17.171 Part (i) of the theorem comes from two facts as 
shown in the following. 

On the one hand, by Lemma 17.201 as soon as z G C is a point in the support 
of Z, we have D(0, \z\) C Supp(Z), where -D(0, \z\) denotes the open disc with 
center and radius \z\. 

On the other hand, the support of Z is unbounded. Indeed, as in (I7.2p . at the 
beginning of the proof of Lemma [7.201 as soon as z G C \ {0} is a point in the 
support of Z, for any t > and for any n G N, m n e~ xt z is in the support of Z. 

For Part (ii), notice that by Lemma [7.25[ for all a G]0, [■> 

ip(t) = 0(t~ a ) as |*| -> oo. (7.6) 

Since EZ ^ 0, by Eq. (15.121) we obtain mEe~ XT = 1, hence mE|e~ AT | = 
mEe -s»(A)T > ! as soon as ^x) ^ 0. Notice that if 3(A) = 0, then A = 1 
by the equation mEe~ XT = 1. So the hypotheses A ^ 1 and EZ ^ imply that 
3(A) ^ and 9fc(A) < 1 (cf. ([5TB]) ). It follows that (JZSJ holds for some a > 1, so 
that the Fourier transform tp of Z is in L 2 . Therefore by the inversion formula 
of Fourier-Plancherel transform, the distribution of Z has a density in L 2 with 
respect to the Lebesgue measure on C. This ends the proof of Theorem 17.171 □ 

Remark 7.26 In fact we have the following more general result. Let A be an 
complex number with a := 3?(A) > 0, r := 3(A) ^ and satisfying the arithmeti- 
cal condition: 

£ Q 



r logm 

and let Z be a nontrivial solution of Eq. A5.11\) (with or without first moment). 
Then the distribution of Z is absolutely continuous with respect to the Lebesgue 
measure on C, and its support is the whole complex plane C. 

To see the conclusions of Remark 17.261 we can argue as follows. In the general 
case where the expectation of Z may not exist, Lemma 17.201 still holds thanks 
to the arithmetical condition. The remaining of the proof is the same, except at 
the end, where 3?(A) > 1 is no more ensured. Nevertheless we have an additional 
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argument by iteration. Iterating the distributional equation ( 15.111) . we obtain for 
n > 1, 

Z — V 1 A A A y(u 1 ...u n ) 

ui...«»e{l,...,m} n 

where A = e~ XT , A u are independent copies of A (indexed by finite sequences of 
integers u), Z^ are independent copies of Z, the two families {A u } and {Z^} 
are also independent of each other; by convention, A Ul . . . A Ul „ Un _ l is taken to be 
1 when n = 1. It is convenient to rewrite this equation in the form 

Z k YY^\ (7.7) 

3=1 

where Z^ are independent copies of Z which are also independent of {Yj}. For 
fixed y = {yj : 1 < J < m n ) with n^L™ Vj ^ 0, by Lemma OS for a e]0, l/3ft(A)[ 
and some constant c > 0, 



Eexp(*(t,^^)» 

3=1 



m 



3=1 

C(y)\t\ 



where C(y) = YYjLi c l%'l a > does not depend on t. Let n > 1 be large enough 
such that m n a > 1. Then X^Li VjZ® is absolute continuous (with respect to the 
Lebesgue measure on C) as its Fourier transform is square integrable on C. This 
implies that for each Borel set B of C with Lebesgue measure 0, we have 

3=1 

It follows from Eq. 17.71 (by conditioning on (Yj)) that P(Z e B) = 0. 

8 Exponential moments and Laplace series 

In this section we consider a solution Z of Eq. (15. lip and we show that its expo- 
nential moments exist in a neighborhood of 0, so that the moment exponential 
generating series of Z defines an analytic function in a neighbourhood of the 
origin. We show that this function satisfies a very simple differential equation. 
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Theorem 8.27 Let A G C be a root of the characteristic polynomial Ii2.3\) with 
3ft (A) > 1/2 and let Z be a solution of Eq. Ii5.11\) . There exist some constants 
C > and e > such that for all t G C with \t\ < e, 

Ee (t,Z) < e «(t)+C|t|» md Ee \tZ\ < 4e |t|+2C|t|»_ (g-Q 

To prove this theorem, we use Mandelbrot's cascades in the complex setting 
(see Barral et al. [2] for independent interest about complex Mandelbrot's cas- 
cades). We still denote A = e~ XT . Then mKA = 1 because A is a root of the 
characteristic polynomial (12 .3p and mE|v4| 2 < 1 because 3ft(A) > 1/2 (see (I5.6P ). 
Let A u ,u G U be independent copies of A, indexed by all finite sequences of 
integers 

u = Ul ...u n G U := |J{l,2,...,m} n 

n>l 

and set Yq = 1, Y\ = mA and for n > 2, 

^ ^AA Ul A UlU2 . . . A Ul _ Un _ 1 . 

As mEA = 1, (Y n ) n is a martingale. This martingale has been studied by many 
authors in the real random variable case, especially in the context of Mandelbrot's 
cascades, see for example [T3] and the references therein. It can be easily seen 
that 

m 

Y n+1 = Aj2Y n ,i (8.2) 

i=l 

where Y n ^ for 1 < i < m are independent of each other and independent of A and 
each has the same distribution as Y n . Therefore for n > 1, Y n is square-integrable 
and 

VarF n+1 = (E\A\ 2 m 2 - 1) + mE\A\ 2 V&rY n , 

where Var X = E (\X — EX| 2 ) denotes the variance of X. Since mE|A| 2 < 1, the 
martingale (Y n ) n is bounded in L 2 , so that the following result holds. 

Lemma 8.28 Let X be a root of the characteristic polynomial k2. 3\) with 3ft(A) > 
1/2. Then, when n — > +oo, 

Y n — > Yoo a.s. and in L 2 , 
where Y^ is a (complex-valued) random variable with variance 



Notice that, passing to the limit in ( 18. 2 j) gives a new proof of the existence 
of a solution Z of Eq. (15.111) with EZ = 1 and finite second moment whenever 
3?(A) > 1/2. From Section [61 we have the unicity of solution of this equation so 
that Theorem 18.271 is proved as soon as it holds for Y^. 

Lemma 8.29 Under the condition of Lemma \ 8.28\ there exist some constants 
C > and e > such that for all t G C with \t\ < e, we have 



Ee {t ' Yoo) < e m+c{t 



(8.3) 



Proof. As in [22] and [16] (where a similar problem for real random variables 
was considered), we use an induction argument. Notice that by Eq. (18. 2p . writing 



<p n (t) :=Ee<*' yn >, teC,n>0, 



we have 



ip n+1 {t)=E<p™{At), 



t g C. 



U) 



We shall prove that there exist some constants C > and e > such that for all 
n > and all £ G C with \t\ < e, we have 



<Pn(t) < e^ +c ^' 



•5) 



Let us prove ( 18. 5 p by induction. The inequality holds clearly for n = since 
Vo(^) = e 3 *^. Assume that it holds for some n > and all t G C with \t\ < e. 
Then writing A = Ai + iA 2 (Aj G R), using |A| < 1 and Eq. ( 18. 4p . we have for 
t = t 1 + it 2 (U G R) with \t\ < e, 



5.6) 



¥>n+i(*) < Eexp{m(Ai*i + A 2 t 2 + C\A\\t{ + t 2 2 ))} 

where g(t 1 ,t 2 ) = Ee^ 1 '^ with 

t 2 ) = (mA 1 - + mA 2 t 2 + C(m|A| 2 - + ^). 

Notice that (7(0,0) = 1. It remains to prove that (0,0) is a local maximum of g. 
Clearly, 



If = Ee" 
dti 

d 2 g 



dh 
dti 



dU 
d 2 g 



TO ft 



i = 1,2 

d 2 h 



^dtj + dt 



dh dh d 2 h 



1,2 



Notice that, a.s. 



(0,0) = (7^-1), ^(0,0)=m4 
Ot\ ot 2 

(0, 0) = 0, ^(0, 0) = 2C(m\A\ 2 - 1), i = 1, 2. 



dtidt 2 dti 
Recall that mEA = 1, so that mMAi = 1 and mEA 2 = 0; hence 

J|(0, 0) = EMj - 1) = 0, J|(0, 0) = E(mi 2 ) = 0, 

so that (0, 0) is a critical point of g. Moreover, 

§^(0,0) = E[(mA 1 -lf + 2C(m\A\ 2 -l)], 
^(0,0) = E \(mA 2 ) 2 + 2C(m\A\ 2 - 1 



<% 2 
dtidt 2 



(0,0) = E(mAi - l)(mA 2 ). 



As E(m|v4| 2 - 1) < (recall that 3?(A) > 1/2), it follows that the Hessian matrix 
at (0, 0) is definite negative for C > large enough which implies that g(0, 0) is 
a local maximum of g. So for e > small enough, g(ti,t 2 ) < #(0,0) = 1 for all 
t = t± + it 2 with \t\ < e. Hence by ( 18. 6ft . for such e and C which do not depend 
on n, ( 18. 5 p holds for n + 1. Therefore, by induction, it holds for all n > 0. 

Letting n — > oo in (18. 5p . we see that inequality ( 18. 3 p holds by Fatou's lemma. 

□ 



Proof of Theorem 18.271 By the unicity of solution of Eq. ( 15. lip . C(Z) = 
£(Yoo). So by Lemma T8.291 there are some constants C > and e > such 
that the first inequality of (18. ip holds. To show the second one, notice that 
|£||9R(Z)| + |t||3(Z)| takes one of the four values ±|t|3?(Z) ± |t|S(Z) (according 
to the signs of 9ft (Z) and $s(Z)), so that a.s. 

e \tZ\ < e |t||»(Z)|+|t||9(Z)| 

< e \t\St(Z)+\t\Q(Z) + e |t|R(^)-|t|9f(Z) + e -\t\St(Z)+\t\SH(Z) + e -|t|«(Z)-|t|9F(^)_ 

Taking expectation in both sides, and noticing that ±\t\$t(z)±\t\s(z)={(±i±i)\t\,z) , we 
see that the second inequality in (18. ip follows from the first one. □ 
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Suppose that Z is any solution of Eq. ( 15. lip under the assumptions of The- 
orem 18.271 The second inequality ( 18. ip shows that the exponential generating 
series of absolute moments of Z has a positive radius of convergence so that the 
formal Laplace series 

p>0 ^' 

defines an analytic function in a neighbourhood of the origin. One can also write 
L(z) = ~Ee zZ when \z\ is sufficiently small. 

Let's come back to the dislocation equations (15. 7p satisfied by the limit vari- 
ables Wi, . . . , W m . These variables admit finite (absolute) moments at any order. 
For any k 6 {1, . . . , m}, let L k be the formal Laplace series defined by 



L k (z) := Yl 



pi 



The dislocation equations (15. 7p imply recursive relations on W k s moments. De- 
veloping these relations with the multinomial formula implies that L k satisfy the 
formal differential system 



A 2 



Vfce{l,...,m-2}, L k (z) + -±zL' k {z) = L k+1 (z) 
A 

L m _iO) + 2 —zL' m _ 1 (z) = {L l {z)) m , 

m — 1 



.7) 



with boundary conditions 

L fc (0) = 1, 1 < k < m- 1, 
L' k (0)=E(W k ) = u 2 (X fc (0)) 



A 2 + k - V 
k-1 

Since W\ satisfies the assumptions of Theorem 18.271 the series L\ has a positive 
radius of convergence as shown above. Therefore, the same holds for all L k be- 
cause of the system (18 . 71) so that the L k define, near the origin, analytic functions 
related by (18~7j) . 

Let p be any complex (m— l)-th root of (— l) m (m— 1)!. For any k G {1, . . . , m}, 
define 

L k (z- X2 ) 
G k (z) := (-l) k p(k - 1)! \ 

z K 



where z~ Xi denotes any determination of the logarithm. For sufficiently large 
\z\, this formula defines an analytic function on a slit plane. Reporting in for- 
mula ( 18. 7p shows that the functions G k satisfy the simple differential system 

V*e{l,...,m-2}, G' k = G k+1 , 

fit rim 

~ ^1 ■ 



In particular, G\ is solution of the differential equation y( m ^ = y m . We sum up 
these results in the following statement. 

Theorem 8.30 Let W\ be the complex-valued limit distribution for the multitype 
branching process of m-ary search trees as defined in Section [57B . Then: 

(i) the Laplace series L 1 (z) = E (e zH/l ) has a positive radius of convergence; 

(ii) for any determination of the logarithm, the function 

z^- P -L, (*-*») , 

z 

is a solution of the differential equation 

Remark 8.31 As can be straightforwardly checked, the function y K (z) := is 
a solution of Eq. (E2j) when the complex number n satisfies K m ~ x = (m — 1)!. 
Nonetheless, G\ is not a function of this form. 

Indeed, since Li(w) = 1 + w + o(w) in a neighbourhood of the origin, G\ 
admits the expansion 

while y K satisfies 

k k k ( 1 

i = 2 +0 \ — 

1 — z z z z \z z 

One concludes by unicity of (complex) power expansions, because A2 7^ 1. 
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